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INTRODUCTION 

Monitoring  biological  aerosols  is  important  and  has  application  to  DoD,  medical  and  health, 
food  processing  and  environmental  monitoring  communities.  An  optical  Ultra  Violet  (UV) 
fluorescence  detector  has  been  developed  by  the  Optical  Sciences  Division  of  NRL  (Code  5615) 
for  this  purpose.  The  detector  consists  of  an  inlet  nozzle  which  directs  sample  particles  to  the 
focus  of  an  elliptical  mirror.  A  red  laser  beam  illuminates  the  flow  at  this  focal  plane  and  the  scatter 
red  light  triggers  a  pulsed  UV  laser.  The  UV  fluorescence  from  the  aerosol  is  then  detected  and 
then  correlated  to  the  aerosol  concentration. 

The  nozzle  employed  in  the  UV  fluorescence  detector  is  similar  to  a  pitot  tube  in  construction. 
An  inner  capillary  tube  carries  the  particulate  flow  while  the  outer  tube  carries  a  sheath  flow.  The 
outer  tube  ends  in  a  tapered  nozzle  and  the  inner  tube  may  either  be  flush  with  the  exit  or  recessed 
inside  the  nozzle.  In  the  present  study,  we  perform  computation  of  flow  exiting  the  nozzle  to 
determine  the  spreading  of  the  flow  at  a  specified  distance  from  the  exit  of  the  nozzle  where  the  red 
laser  beam  is  directed.  A  series  of  computations  are  performed  varying  the  distance  by  which  the 
inner  tube  is  recessed  and  the  flow  rates  through  the  nozzle  and  study  the  effect  on  the  spreading 
rate  of  the  exiting  flow. 


INCOMPRESSIBLE  FLOW  SOLVER 

FEFLO  is  a  finite  element-based  incompressible  flow  solver  based  on  simple,  low-order 
elements.  The  simple  elements  enable  the  flow  solver  to  be  as  fast  as  possible  reducing  the 
overhead  in  building  element  matrices,  residual  vectors  etc.  The  governing  equations  are  written  in 
Arbitrary  Lagrangian  Eulerian  form  which  enables  simulation  of  flow  with  moving  bodies.  For 
high  Reynolds  number  flow  cases,  the  mesh  requirement  is  met  by  employing  arbitrary  semi- 
structured  grids  close  to  wetted  surfaces  and  wakes.  The  full  details  of  the  flow  solver,  the  rigid 
body  motion  and  adaptive  remeshing  are  given  by  Ramamurti  et  al.1  and  are  summarized  next. 

The  governing  equations  employed  are  the  incompressible  Navier-Stokes  equations  in 
Arbitrary  Lagrangian  Eulerian  (ALE)  formulation.  They  are  written  as 
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Here  p  denotes  the  pressure,  va  =v-w,  the  advective  velocity  vector  (flow  velocity  v minus 

mesh  velocity  w),  and  the  material  derivative  is  with  respect  to  the  mesh  velocity  w.  Both  the 
pressure  p  and  the  stress  tensor  a  have  been  normalized  by  the  (constant)  density  p,  and  are 
discretized  in  time  using  an  implicit  time  stepping  procedure.  It  is  important  for  the  flow  solver  to 
be  able  to  capture  the  unsteadiness  of  a  flow  field.  The  present  flow  solver  is  built  as  time-accurate 
from  the  onset,  allowing  local  timestepping  as  an  option.  The  resulting  expressions  are 
subsequently  discretized  in  space  using  a  Galerkin  procedure  with  linear  tetrahedral  elements.  In 
order  to  be  as  fast  as  possible,  the  overhead  in  building  element  matrices,  residual  vectors,  etc. 
should  be  kept  to  a  minimum.  This  requirement  is  met  by  employing  simple,  low-order  elements 
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that  have  all  the  variables  (u,v,w,p)  at  the  same  location.  The  resulting  matrix  systems  are  solved 
iteratively  using  a  preconditioned  gradient  algorithm  (PCG).  The  flow  solver  has  been  successfully 
evaluated  for  both  2-D  and  3-D,  laminar  and  turbulent  flow  problems  by  Ramamurti  et  al.2  i 


COMPUTATIONAL  RESULTS 

The  nozzle  employed  in  the  UV  fluorescence  detector  is  similar  to  a  pitot  tube  in  construction 
and  is  shown  in  Fig.  1 .  An  inner  capillary  tube  carries  the  particulate  flow  while  the  outer  tube 
carries  a  sheath  flow.  The  outer  tube  ends  in  a  tapered  nozzle  and  the  inner  tube  is  flush  with  the 
exit  of  the  nozzle.  Figure  2  shows  the  variation  of  the  pressure  drop  across  a  nozzle  of  similar 
construction  for  various  flow  rates.  These  results  were  obtained  experimentally  by  J.  Eversole  of 
Code  5615.  Figures  3a-d  show  the  results  of  the  computations  for  a  total  flow  rate,  Qtotal,  of 
12Lpm  with  the  inner  particulate  flow,  Qinner,  of  3Lpm  and  the  outer  sheath  flow,  Qouter,  of  9Lpm. 
The  computational  grid  consists  of  537K  points  and  2.9M  tetrahedral  elements.  The  surface 
triangulation.  Fig.  3a,  shows  the  grid  clustering  near  the  exit  of  the  nozzle.  From  Fig.  3b,  the 
pressure  drop  can  be  seen  to  be  linear  inside  the  nozzle.  The  computational  unit  of  pressure  for  the 

incompressible  flow  is  m2/s2  (=p/p),  where  p  is  the  density  of  air  (Kg/m3).  The  computed 
pressure  drop  for  this  flow  is  47.5torr  and  the  experimental  value  from  Fig.  2  is  approximately 
45.0torr.  Figure  3c  and  3d  show  the  variation  of  velocity  near  the  exit  of  the  nozzle.  It  can  be 
seen  that  the  inner  and  outer  flows  are  nearly  parallel  with  minor  instabilities  present  near  the  edges 
of  the  outer  flow. 

The  nozzle  geometry  was  next  updated  and  is  shown  in  Fig.  4.  The  inner  capillary  tube  in  this 
configuration  has  a  taper  and  the  inner  diameter  of  the  outer  tube  is  also  reduced  to  4.0mm.  The 
experimentally  observed  pressure  drop  for  this  nozzle  is  shown  in  Fig.  5.  Computations  were 
performed  for  three  different  flow  rates,  4Lpm,  8Lpm  and  12Lpm.  The  inner  and  outer  flow  rates 
are  prescribed  for  these  computations.  For  a  QlotaJ  =  4Lpm,  the  prescribed  flow  rates  were  Qinner  = 
1 .5Lpm  and  Qouter  =  2.5Lpm;  for  Qtotal  =  8Lpm,  the  prescribed  flow  rates  were  Qinner  =  2.5Lpm  and 
Qouter =  5.5Lpm;  for  Qtotal  =  12Lpm,  Qinner  =  3.0Lpm  and  Qouter  =  9.0Lpm.  The  computed  pressure 
drops  for  the  inner  and  outer  flows  are  shown  in  Fig.  6.  The  average  pressure  drop  across  the 
entire  nozzle  is  computed  as  an  area  weighted  average  of  the  pressure  drops  of  the  inner  and  outer 
flows.  The  pressure  drop  for  the  inner  flow  increases  from  12.3torr  to  47.4torr  as  the  inner  flow 

rate  Q^ner  *s  increased  from  1.5Lpm  to  3.0Lpm.  For  the  outer  flow,  the  Ap  drops  from  38.65torr 
to  0.7torr  as  Qouter  is  increased  from  2.5Lpm  to  5.5Lpm.  At  Qtotal  =  12Lpm,  the  computational 
pressure  drop  is  47.4torr  compared  to  45.0  in  the  experiments.  The  prescribed  computational  flow 
rates  are  Qinner  =  3Lpm  and  Qouter  =  9Lpm,  whereas  the  experimentally  observed  flow  rates  are 
2.5Lpm  and  8.0Lpm,  approximately.  It  should  be  noted  that  these  flow  rates  do  not  add  up  to  the 
total  flow  rate  of  12.0Lpm.  For  a  total  flow  rate  of  8Lpm,  the  computational  pressure  drop  is 
38.8torr  compared  to  23.0torr  observed  in  the  experiments.  The  prescribed  Qinnei.  and  Qouter  for  this 
case  are  2.5Lpm  and  5.5Lpm,  respectively.  Again,  there  is  a  discrepancy  in  the  experimentally 
observed  inner  and  outer  flow  rates  and  the  total  flow  rate  by  approximately  0.75Lpm.  At  Qtotal  = 
4Lpm,  the  computational  pressure  drop  is  12.3torr  compared  to  8.0torr  observed  in  the 
experiments.  These  differences  may  be  due  the  experimental  errors  in  measuring  the  flow  rates 
and  the  experimental  error  is  not  known.  The  prescribed  Qinner  and  Qouter  for  this  case  are  1.5Lpm 
and  2.5Lpm,  respectively. 

Next,  the  position  of  the  exit  of  the  inner  capillary  tube  is  recessed  by  5.0mm  and  10.0mm  the 
computations  are  performed  for  these  configurations  for  the  three  flow  rates.  Figures  7a-c  show 
the  variation  of  the  centerline  pressure  for  the  three  different  nozzle  configurations.  It  can  be  seen 
that  as  the  flow  rate  is  increased,  the  pressure  drop  across  the  nozzle  increases,  as  would  be 
expected.  The  effect  of  the  capillary  recession  on  the  centerline  pressure  distribution  is  shown  in 
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Figs.  8a-c.  For  a  flow  rate  of  4Lpm,  Fig.  8a,  the  pressure  drop  for  the  case  where  the  capillary  is 
recessed  by  5mm  exceeds  that  of  the  flush  and  the  10mm  cases.  This  trend  is  reversed  as  the  flow 
rate  is  increased  to  8Lpm  and  12Lpm.  For  all  the  three  flow  rates  studied,  the  pressure  drop  for 
the  10mm  case  is  higher  than  the  flush  capillary.  In  order  to  understand  the  reason  for  the  trend 
reversal  for  the  5mm  case,  the  details  of  the  flow  for  these  cases  are  examined.  Figures  9  and  10 
show  the  variation  of  pressure  drop  with  flow  rate  when  the  inner  capillary  tube  is  recessed  by 
5.0mm  and  10.0mm,  respectively.  From  Figure  9,  it  can  be  seen  that  the  pressure  drop  across  the 
inner  capillary  increases  as  the  flow  rate  is  increased  from  1.5Lpm  to  2.5Lpm  and  remains  almost 
same  for  a  flow  rate  of  5.5Lpm.  The  pressure  drop  for  the  outer  flow  behaves  similar  to  that  of  the 

flush  capillary  case.  Figure  10  shows  that  the  Ap  increases  almost  linearly  with  flow  rate. 

Figure  1 1  shows  the  surface  definition  and  the  triangulation  of  the  nozzle  geometry  with  flush 
capillary.  The  computational  mesh  consists  of  242K  points  and  1.2M  tetrahedral  elements.  Figure 
12  shows  the  results  of  the  computation  for  this  configuration  at  a  flow  rate  of  12Lpm.  The 
pressure  drop  inside  the  nozzle  is  almost  linear,  Fig.  12a,  and  is  almost  constant  after  the  exit  of 
the  nozzle.  The  flow  exiting  the  nozzle  is  parallel  as  can  be  seen  from  the  velocity  contours, 
velocity  vectors  and  particle  traces.  Figs.  12b-d.  Figure  13a  shows  the  particle  traces  for  the 
nozzle  with  flush  capillary  at  4Lpm,  and  the  flow  exiting  is  nearly  parallel,  and  the  outer  sheath 
flow  exhibits  a  divergence,  which  is  not  of  primary  interest  of  the  study.  Figure  14  shows  the 
geometry  of  the  nozzle  with  the  capillary  tube  recessed  from  the  exit  of  the  nozzle.  Figures  15a 
and  b  show  the  pressure  inside  the  nozzle  and  the  velocity  vectors  near  the  exit,  for  a  nozzle  with 
the  capillary  recessed  by  5.0mm.  It  is  clear  that  at  the  edge  of  the  inner  particulate  flow,  shear 
layer  instability  is  present.  This  is  also  evident  in  the  particle  traces  shown  in  Figs.  16a-c,  and  is 
present  for  all  the  three  flow  rates  studied.  Figure  17  shows  the  contours  of  the  magnitude  of 
velocity  outside  the  capillary  tube  and  the  nozzle  exit.  The  capillary  tube  is  recessed  by  10.0mm  in 
this  case,  and  the  flow  exiting  it  is  nearly  parallel.  Figures  1 8a-c  show  the  particle  traces  for  this 
configuration  at  three  different  flow  rates.  In  all  the  cases,  the  particulate  flow  is  nearly  parallel, 
and  the  sheath  flow  shows  a  small  instability  at  the  outer  edge  of  the  shear  layer.  Figures  19-21 
show  the  particle  traces  for  the  three  different  nozzle  configurations  and  for  the  minimum  and 
maximum  flow  rates  considered  in  this  study.  The  particles  traces  shown  are  released  at  0.2mm, 
0.4mm  and  0.6mm  from  the  centerline  of  the  nozzle.  For  this  nozzle,  the  inner  radius  of  the 
capillary  is  0.375mm  and  that  of  the  nozzle  is  2.0mm,  Fig.  4.  For  the  nozzle  with  the  flush 
capillary  tube,  Figs.  19a  and  b,  it  can  be  seen  that  the  inner  flow  is  convergent  in  the  vicinity  of  the 
nozzle  exit.  Particles  released  at  0.6mm  stray  away  from  the  centerline  in  the  vicinity  of  the  exit 
and  are  drawn  to  the  core  flow  downstream.  As  the  flow  rate  is  increased,  the  tendency  for  the 
outer  flow  to  spread  is  reduced.  For  the  case  of  the  nozzle  with  the  capillary  recessed  by  5.0mm, 
Figs.  20a  and  b,  the  inner  the  flow  exhibits  shear  layer  instability.  In  this  case  also,  the  particles 
released  at  0.6mm  tend  to  spread  more  at  lower  flow  rates  and  are  drawn  to  the  inner  flow  as  the 
flow  rate  is  increased.  For  the  nozzle  configuration  with  the  capillary  recessed  by  10.0mm,  Figs. 
21a  and  b,  the  inner  flow  remains  parallel  and  the  particles  released  at  0.6mm  moved  closer  to  the 
inner  flow  as  the  flow  rate  is  increased.  The  particle  traces  in  the  region  illuminated  by  the  red 
laser  beam  is  shown  in  Figs.  22a  and  b.  At  lower  flow  rate  of  Qtota]  =  4Lpm,  Fig.  22a,  the  particle 
released  at  0.4mm  shows  a  small  instability  for  the  flush  capillary  and  the  5.0mm  recessed  case. 
As  the  capillary  is  recessed  even  more,  this  instability  disappears,  but  the  spreading  rate  of  the 
inner  flow  is  increased.  At  a  higher  flow  rate  of  12Lpm,  Fig.  22b,  the  instability  is  not  evident. 
As  the  flow  rate  is  increased,  the  inner  flow  becomes  more  convergent  for  the  nozzle  with  the 
flush.  For  the  nozzle  with  the  capillary  recessed  by  10.0mm,  the  spreading  rate  seems  to  be 
independent  of  the  flow  rate. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  incompressible  flow  solver  FEFLO  was  used  to  compute  the  flow  through  the  nozzle  of  a 
UV  Fluorescence  detector.  Computed  pressure  drop  for  the  nozzle  is  compared  with  experiments 
and  is  in  good  agreement.  Three  different  nozzle  configurations  were  studied  by  varying  the 
distance  by  which  the  inner  capillary  tube  is  recessed.  For  each  of  these  configurations,  the  flow 
rate  through  the  nozzle  was  varied.  The  spreading  rate  of  the  flow  exiting  the  nozzle  was 
computed.  For  the  nozzle  with  the  flush  capillary  tube,  the  inner  flow  remains  convergent  and  as 
the  flow  rate  is  increased  for  this  configuration,  it  also  becomes  stable.  For  the  nozzle  with  the 
capillary  recessed  by  5.0mm,  the  inner  and  outer  flows  exhibit  shear  layer  instability,  for  all  the 
flow  rates  considered,  and  for  this  reason  may  not  be  desirable.  The  nozzle  with  the  capillary 
recessed  by  10.0mm  shows  the  least  instability  of  the  inner  flow,  although  the  spreading  rate  is 
higher  than  the  other  two  cases.  The  spreading  rate  is  also  nearly  independent  of  the  flow  rate,  and 
may  be  the  most  desirable  of  all  the  three  cases  considered  in  this  study. 
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b.  Pressure 
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d.  Velocity  vectors  near  the  exit 
Figure  3.  Flow  through  the  nozzle,  Qtotal  =  12Lpm 
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Ap  (torr.) 

Figure  6.  Computed  pressure  drop 
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a.  Flush  capillary 

Figure  7.  Effect  of  flow  rate  on  centerline  pressure 
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a.  Surface  definition 


b.  Surface  mesh 


Figure  1 1 .  Nozzle  geometry 
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d.  Particle  traces 

Figure  12.  Flow  through  the  nozzle  with  flush  capillary,  Qtoul  =  12Lpm 
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a.  Pressure 

Figure  15.  Flow  through  the  nozzle  with  capillary  moved  in  5.0mm,  Q[otal  =  4Lpm 
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b.  Velocity  vectors 

Figure  15.  Flow  through  the  nozzle  with  capillary  moved  in  5.0mm,  Qtotal  =  4Lpm 


a-  Qtota,  =  4Lpm 

Figure  16.  Particle  traces  for  nozzle  with  capillary  moved  in  5.0mm 
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c.  Q,ot,i  =  12Lpm 

Figure  16.  Particle  traces  for  nozzle  with  capillary  moved  in  5.0mm 
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Figure  17.  Flow  through  nozzle  with  capillary  moved  in  10.0mm, 
Qtoiai  =  4Lpm,  absolute  velocity  contours 


a.  Qlotal  =  4Lpm 

Figure  18.  Particle  traces  for  nozzle  with  capillary  moved  in  10.0mm 
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b.  QtotaI  =  8Lpm 


c.  Qtota1  =  12Lpm 

Figure  18.  Particle  traces  for  nozzle  with  capillary  moved  in  10.0mm 
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Figure  22.  Effect  of  capillary  position  on  particle  traces 


